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ABSTRACT 

We study the clustering properties of if-selected galaxies at 2 < z < 3.5 using deep multiwavelength 
imaging in three fields from the MUSYC survey. These are the first measurements to probe the 
spatial correlation function of if-selected galaxies in this redshift range on large scales, allowing for 
robust conclusions about the dark matter halos that host these galaxies. iC-selected galaxies with 
K < 21 have a correlation length rg ~ 6ft.^^Mpc, larger than typical values found for optically- 
selected galaxies. The correlation length does not depend on K-hand magnitude in our sample, but it 
does increase strongly with color; the J — K > 2.3 distant red galaxies (DRGs) have tq ~ ll/i~^Mpc. 
Furthermore, contrary to findings for optically-selected galaxies i^T-selected galaxies that are faint in 
the _R-band cluster more strongly than brighter galaxies. These results suggest that a color-density 
relation was in place at z > 2; it will be interesting to see whether this relation is driven by galaxies 
with old stellar populations or by dusty star forming galaxies. Irrespective of the cause, our results 
indicate that if-bright blue galaxies and ii'-bright red galaxies are fundamentally different, as they 
have different clustering properties. Using a simple model of one galaxy per halo, we infer halo masses 
'--^ 5 X lO^^M© for K < 21 galaxies and ~ 2 x IO^^Mq for DRGs. A comparison of the observed space 
density of DRGs to the density of their host halos suggests large halo occupation numbers; however, 
this result is at odds with the lack of a strong small-scale excess in the angular correlation function. 
Using the predicted evolution of halo mass to investigate relationships between galaxy populations at 
different redshifts, we find that the z = descendants of the galaxies considered here reside primarily 
in groups and clusters. 

Subject headings: galaxies: evolution — galaxies: formation — galaxies: high-redshift — cosmology: 
large-scale structure of the universe — infrared: galaxies 



1. INTRODUCTION 

Optical surveys of the high-redshift universe have been 
very successful in finding relatively unobscured star- 
forming galaxies, primarily via the [/-dropout technique. 
These z ~ 3 Lyman break galaxies (LBGs) typically 
have stellar masses ^ IO^'^^Mq, star formation rates of 
10 — lOOM0yr~^, and are thought to dominate the star 
formation density at that epoch (Steidel et al. 2003; 
Shapley et al. 2001; Reddy et al. 2005). However, it 
is becoming increasingly clear that substantial numbers 
of galaxies exist at these redshifts that have little rest- 
frame UV luminosity and are thus underrepresented in 
optical surveys. Such galaxies may be detected in the 
near-infrared (NIR), which samples the rest-frame opti- 
cal. 

One criterion used to select galaxies in the NIR is 
J - K > 2.3 (Franx et al. 2003; van Dokkum et al. 
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2003). These distant red galaxies (DRGs) typically have 
high star formation rates, > lOOMoyr"^, and dust ob- 
scuration Ay > 1, but must also have significant pop- 
ulations of evolved stars in order to explain their colors 
and spectra (Forster Schreiber et al. 2004; Papovich et 
al. 2006; Kriek et al. 2006a). Some DRGs show little or 
no evidence of active star formation (Labbe et al. 2005; 
Kriek et al. 2006a,b; Reddy et al. 2006). DRGs must 
in general be very massive to account for their signif- 
icant X-brightness; stellar population synthesis models 
imply masses ^ lO^^M©. Indeed, 95% of galaxies 
with Af, > 1O"M0 at 2 < z < 3 have K < 21.3 (van 
Dokkum et al. 2006). Conversely, the median galaxy in 
this mass range has R ~ 25.9, fainter than the limits 
typically reached by optical surveys. 

The relationship between i^T-selected samples and 
optically-selected samples, not to mention present-day 
galaxies, remains unclear. While typical optically- 
selected galaxies have different properties than K- 
selected galaxies, the i^T-bright subsample of optically- 
selected galaxies have stellar masses, star formation 
rates, and metallicities that are in approximate agree- 
ment with X-selected galaxies (Shapley et al. 2004; 
Reddy et al. 2005). Nevertheless, it is not clear whether 
the differences between these galaxies are transient (e.g. 
dust geometry, starbursts, mergers) or fundamental (e.g. 
age of the underlying old stellar populations, mass, envi- 
ronment) . An understanding of the nature of the differ- 
ences between populations is essential to place them in 
evolutionary scenarios. 
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One way to investigate differences between galaxy pop- 
ulations is to measure their clustering properties. As 
clustering measurements provide information that is in- 
dependent of photometric properties, they can be used 
to distinguish between transient and fundamental differ- 
ences between galaxy populations. In the halo model of 
galaxy formation, the large-scale distribution of galaxies 
is determined by the distribution of dark matter halos. 
The correlation function of galaxies can therefore be asso- 
ciated with the correlation function of the halos in which 
they reside. Halo clustering, in turn, is a strong function 
of halo mass (Mo & White 1996), providing a means to 
study the relationship between galaxy properties and the 
mass of the dark matter halos. 

Several studies have measured the dependence of clus- 
tering strength of high-redshift galaxies on color. Daddi 
et al. (2003) use the ultradecp imaging of the 4.5arcmin^ 
FIRES HDF-S field (Labbe et al. 2003) to study the 
clustering characteristics of iiT-selected galaxies at 2 < 
z < 4. Their most striking finding is that the corre- 
lation length increases strongly with J — K color, with 
the reddest galaxies in their sample having correlation 
lengths To = 10 — 15/i~^Mpc, comparable to the most 
luminous red galaxies in the local universe. The mass 
of dark matter halos with similar correlation lengths is 
> 1O"'^^M0, yet the galaxy number density is ~ 100 times 
larger than that expected for such massive dark matter 
halos, implying that many galaxies must share the same 
halo. More recently, Grazian et al. (2006) measured a 
ro = 13.4t3;2/i"^Mpc for z > 2 DRGs using the larger 
135arcmin2 GOODS-CDFS field, also indicating that red 
galaxies are located in very massive halos. 

The interpretation of these clustering measurements 
is complicated by the fact that, in order to derive in- 
formation about the dark matter halos, the correlation 
function must be measured on large scales. The correla- 
tion function has a contribution from galaxies that share 
halos (hereafter the "1-halo" term) and galaxies in sep- 
arate halos (the "2-halo" term). The shapes of these 
two contributions are shown with impressive detail in 
the correlation functions of large samples of z ~ 4 LBGs 
presented by Lee et al. (2006) and Ouchi et al. (2005). 
In order to derive meaningful constraints on large-scale 
clustering properties, and thus on the host dark matter 
halos, it is important that both of these terms are taken 
into consideration. If the correlation function is parame- 
terized as a simple power-law then it should be measured 
on scales where the 2-halo term dominates the clustering 
signal. In practice, a firm lower limit to the radial range 
in which the angular correlation function should be fit- 
ted is the halo virial radius. At z = 3, the virial radius 
r2oo of a lO"'^''Af0 halo corresponds to 22" (e.g. Mo & 
White 2002). The majority of the clustering signal from 
Daddi et al. (2003) and Grazian et al. (2006) is on scales 
< 30", which may lead to gross overestimates of the 
large-scale correlation length and the mass of the host 
dark matter halos. In particular, Zheng (2004) shows 
that the measurements of Daddi et al. (2003) are con- 
sistent with models in which the large-scale correlation 
length is as low as ro ~ 5h~^ and the typical halo mass 
is ~ IO^^Mq. This suggests that DRGs and LBGs may 
occupy similar halos, but that DRGs have higher occu- 
pation numbers. 



The goal of this work is to study the clustering char- 
acteristics of a ii'-selected population of galaxies at 2 < 
z < 3.5. The increased field-of-view of our imaging allows 
for an improved determination of the clustering strength 
of i^-selected galaxies at angular separations sufScicntly 
large to investigate the large-scale distribution of galax- 
ies, and thus to provide more meaningful estimates of 
the masses of the halos in which they reside. Secon- 
darily, we wish to analyze the clustering results using 
models of halo clustering, and to use these models to 
shed light on evolutionary scenarios for z > 2 galax- 
ies. These types of analyses have previously been per- 
formed for optically-selected samples (e.g. Moustakas & 
Somerville 2002; Ouchi et al. 2004; Adclberger et al. 
2005a). Throughout, we use the cosmological parame- 
ters 51,„ = 0.3, Qa = 0.7, Ho = IQhm km s^^ Mpc"\ 
and (Tg = 0.9. Results are given using hja = 1, except for 
correlation lengths which are scaled to units of /i = 0.7 in 
order to facilitate comparison to previous studies. Opti- 
cal magnitudes are given in AB and NIR magnitudes are 
given on the Vega system. 

2. DATA 

The Multiwavelength Survey by Yale-Chile (MUSYC) 

consists of optical and NIR imaging of four independent 
30' X 30' fields plus spectroscopic follow-up (Gawiser et 
al. 2006, R. Quadri et al. 2006, in preparation).^ Deeper 
JHK imaging was obtained over 10' x 10' sub-fields with 
the ISPI camera at the CTIO Blanco 4m telescope. The 
present analysis is restricted to three of these deep fields 
(the adjacent HDFSl and HDFS2, and SDSS 1030). The 
deep JHK data will be described in detail elsewhere 
(R. Quadri et al. 2006, in preparation). The total 5(7 
point source limiting depths are J 23.0, H 21.8, 
and K ^ 21.3. The optical UBVRIz data are described 
in Gawiser et al. (2006). 

In this study w(! use spectroscopic redshifts where pos- 
sible, but must rely primarily on photometric redshifts. 
Photometric redshifts were determined using the meth- 
ods of Rudnick et al. (2001, 2003). Briefly, non-negative 
linear combinations of galaxy templates are fit to the 
observed spectral energy distributions. The templates 
include the four empirical templates of Coleman, Wu, & 
Weedman (1980), as well as the two empirical starburst 
templates of Kinney et al. (1996), all of which have been 
extended into the UV and NIR using models. As the 
empirical templates are derived from low-redshift sam- 
ples, we find that they do not adequately describe all 
z > 2 galaxies. For this reason we added lOMyr and 
IGyr old single stellar population templates generated 
with the Bruzual & Chariot (2003) models. The redshift 
probability distribution for each galaxy is calculated us- 
ing Monte Carlo simulations in which the observed fluxes 
are varied within the photometric uncertainties. 

A comparison of the photometric redshifts to spectro- 
scopic redshifts drawn from the literature and from our 
own observations yields a mean Az/(l + z) = 0.12 for 
z > 1.5, corresponding to Az ~ 0.4 at z ~ 2.5. The 
dashed curves in Fig. 1 show the redshift distribution for 
all MUSYC galaxies with 2 < z < 3.5, and for DRGs 
in the same redshift range. The solid curves show the 
distributions that are inferred by summing the redshift 

8 www.astro.yale.edu/MUSYC 
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Fig. 1. — The inferred redshift distribution of galaxies selected 
with 2 < 2:phot < 3.5. The upper two curves are for the full sample, 
and the lower two curves are for galaxies that meet the J — K > 2.3 
criterion for DRGs. Dashed curves indicate the distribution of 
Zphot values, smoothed with a A(2) = 0.4 boxcar average. Solid 
curves indicate the distributions derived by summing the redshift 
probability distributions for each galaxy. The normalization is ar- 
bitrary. 

probability distributions. All distributions have been 

smoothed with a Az = 0.4 boxcar to limit spikes that, 
given our uncertainties, may not be real. 

We restrict the sample to galaxies with 2 < Zphot < 3.5 
and K < 21, except where noted. Fig. 2 shows the po- 
sitions of the DRGs in the three MUSYC fields. Also 
shown are galaxies with stellar mass > IO^^Mq, 
where is determined with stellar population synthesis 
models (§3.4.4). 

3. THE TWO-POINT CORRELATION FUNCTION 
3.1. Method 

The two-point correlation function can be measured by 
counting the number of unique galaxy pairs as a function 
of separation, and comparing the resultant distribution 
to that of a catalog of random points with the same num- 
ber density and subject to the same observing geometry. 
Several estimators for the angular two-point correlation 
function are available, but the estimator introduced by 
Landy & Szalay (1993) is emerging as the de facto stan- 
dard for high-redshift studies. It has been shown to min- 
imize the variance and biases associated with other esti- 
mators (Landy & Szalay 1993; Hamilton 1993; Kerscher 
et al. 2000). The observed amplitude of the two point 
correlation function is thus 



WobsiO) = 



DD{0) - 2DR{0) + RR{e) 
RR{0) ' 



(1) 



where DD{6) is the number of data-data pairs with an- 
gular separation in the interval {6 — + AO/2). 
DR{0) is the number of data-random pairs, and RR{0) 
is the number of random-random pairs, in the same angu- 
lar interval. We use AO = 20". In order to better sample 
the observing geometry, and to decrease the uncertainty 
in DR{6) and RR{9), we use 100 times more random 



points than data points. This requires normalizing the 
DR and RR terms such that DR{e) = J2e ^^i^) = 
EoDDiO). 

The angular correlation function can be approximated 
as a power law 

w{9) = A^0-^. (2) 

However, as the (suitably normalized) number of random 
pairs is equal to the number of data pairs, and since the 
two-point correlation function is the excess probability 
of finding a data pair versus a random pair, it is clear 
that Wobs{(^) cannot be positive for all 9. In particular. 



(3) 



This integral constraint requires that Wobs{9) fall below 
the intrinsic w{0) (Groth & Peebles 1977). The size of 
this bias increases with the clustering strength and de- 
creases with field size; in practice, it is a significant effect 
and a correction must be made. The integral constraint 
correction is approximately constant and equal to the 
fractional variance of galaxy counts in a field. 



IC . 



< Ngal > 



(4) 



where the first term on the right is the Poisson vari- 
ance and the second accounts for the additional variance 
caused by clustering 

"^1 = -^ J J w{di2)dnidQ2 (5) 

(Peebles 1980, §45). Although the clustering term dom- 
inates the integral constraint, the Poisson term is non- 
negligible for the small sample sizes considered here. Fol- 
lowing Infante (1994) and Roche et al. (1999) the clus- 
tering term cr^ can be estimated numerically using 



(6) 



The quantity cr^/A^ is estimated directly from the ran- 
dom catalog for an assumed value p. The amplitude 
of the angular correlation function is related to the ob- 
servations through the fitting function 



wobs{e) = A^e-i^ - IC. 



(7) 



We estimate Ayj iteratively using eqs. 6, 4, and 7. The 
final result is robust against differences in the initial esti- 
mate of Ay^, and convergence only takes a few iterations. 

In the weak clustering regime the uncertainty in the 
Landy & Szalay estimator can be estimated by assuming 
that DD{6) has Poisson variance (Landy & Szalay 1993); 
in this case 

1 + w{0) 



6Wobs{0) 



^/DD{0)' 



(8) 



If the angular correlation function is a power law, the 
spatial correlation function will also be a power law 



(9) 



where tq is the spatial correlation length and j = P + 1. 
The angular correlation function can be used to obtain 
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Fig. 2.— The position of < 21 galaxies at 2 < Zphot < 3.5 in the deep MUSYC fields. The field sizes are ~ 20' X 10' and ~ 10' X 10' 
for HDFSl/2 and 1030, respectively. The large black circles represent galaxies with stellar mass Af* > lO^^M© and the small black circles 
represent less massive galaxies. Galaxies that meet the J — K > 2.3 criterion for distant red galaxies are marked with an open red circle. 



the spatial correlation function by inverting the Limber 
projection 



H^r^^ J F{z)rl'-<{z)N^{z)E{z)dz 
{c/Ho){jN{z)dzy ' 



(10) 



where fc(z) is the comoving radial distance, N(z) is the 
redshift distribution. 



= r(i/2) 



r[(7-i)/2] 
r(7/2) ' 



and 



(11) 



(12) 



(e.g. Maghocchetti & Maddox 1999). 

The function F(z) describes the evolution of clustering 
with redshift, ^(r, z) = ^(r, 0)F{z). The evolution has of- 
ten been modelled as F{z) = (1 + z)^(^~'''+'^', where the 
parameter e is typically specified using e = 7 — 3 for con- 
stant clustering in comoving units, e = for 'stable clus- 
tering', or e = 7 — 1 for 'linear growth' (e.g. Moscardini 
et al. 1998; Overzier et al. 2003). We assume constant 
clustering in comoving units over 2 < 2 < 3.5; this sets 
F{z) = 1. Different values of e, where the correlation 
length is then determined at the median redshift of the 
observed sample, yield similar results. 

3.2. Measurement Strategies 

In what follows we restrict the analysis to galaxies with 
photometric redshift 2 < Zphot < 3.5. Reducing the red- 
shift range produces comparable correlation lengths, but 
with larger uncertainties. Additionally, as we largely rely 
on photometric redshifts, we cannot be confident in our 
ability to divide the sample too finely in redshift space. 

It is common practice in the literature to assume 
7 = 1.8 if the data are not sufficient to make indepen- 
dent measurements of both the slope and the amplitude 
of the correlation function. Recent studies have found 
that 7 ~ 1.6 may be more appropriate for LBGs (Adel- 
berger et al. 2005a; Lee et al. 2006). Direct comparisons 
of the correlation length from different studies can be 
problematic unless the same /3 was used; for this reason, 
the results summarized below use 7 = 1.8 but Tables 1 
and 2 also gives the tq values corresponding to 7 = 1.6. 

In placing the random objects on the image, we mask 
out regions where galaxies could not be detected, e.g. 



in the vicinity of bright stars. This procedure makes a 
negligible change in the resultant correlation functions. 

We measure Wobs{0) (eq. 1) in linearly-spaced 20" bins 
for the purpose of computing the fits. Spacing the 
bins at equal logarithmic intervals gives similar results 
for most samples considered here. We present the results 
of power law fits over the range Q" < 9 < 200" in Table 1, 
and 40" <6 < 200" in Table 2. Most of the discussion 
and analysis in this paper is based on the latter fits. The 
200" upper limit minimizes the effects of low level biases 
(such as errors in the flat-fielding and the integral con- 
straint) and edge effects. The 40" lower limit is set so 
that the fit is not strongly affected by 1-halo term of the 
correlation function (§4). At z = 2, the minimum red- 
shift of our sample, 40" subtends 0.7/i~^Mpc in comoving 
units, which corresponds to the virial radius r2oo of an 

3 X 1O^^M0 halo. This is roughly the mass of halos 
that host the most clustered galaxies in our sample (§4), 
and is larger than the scales (< 10"; < 0.25/i~^Mpc) at 
which z ~ 4 LBGs show significant contributions from 
the 1-halo term (Ouchi et al. 2005; Lee et al. 2006). For 
reference, Adelberger et al. (2005a) uses a lower limit 
of 60", but our reduced signal-to- noise does not allow 
for such a conservative limit. In contrast, Grazian et al. 
(2006) fit w{e) over I" < 6 < 100". Neither Grazian et 
al. (2006) nor Daddi et al. (2003) significantly constrain 
w{9) beyond 40" (see their Figs. 9 and 8, respectively). 

We note that, for most samples considered here, the 
Q" < 9 < 200" fits result in correlation lengths that 
are larger than the 40" < 9 < 200" correlation lengths 
by '-^ 1 — 1.5a. For some subsamples the 0" < 9 < 
200" correlation lengths are actually smaller, although 
the difference is always < la. 

While performing the fit at large scales reduces the ef- 
fect of the 1-halo term on the correlation length, there 
is a second-order effect of the halo occupation distribu- 
tion that we do not take into account. A fully consis- 
tent treatment would require counting only one galaxy 
per halo, to avoid counting the same halo more than 
once. As we have no robust method to detect galaxies 
that share halos, halos that host multiple galaxies will 
be counted multiple times when measuring w{9). Since 
only the most massive halos are likely to host multiple 
galaxies, these halos effectively receive more weight. Our 
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data are not sufficient to address these second-order ef- 
fects, and we note that simply rejecting all galaxies that 
have close neighbors would introduce other biases in our 
measurements. 

3.3. Sources of uncertainty 
3.3.1. Redshift distribution 

The shape of the redshift distribution N{z) in eq. 10 
will affect the deprojection of the angular two point cor- 
relation function, contributing to the uncertainty in ro- 
One strategy for dealing with the redshift uncertainties 
is to smooth the Xphot distribution by the typical un- 
certainty, Az ^ 0.4 (§2). However, it is likely that the 
redshift uncertainties differ for different galaxies in our 
sample. This may affect the observed relationships be- 
tween galaxy properties and clustering strength; for in- 
stance, if faint galaxies have a larger photometric redshift 
uncertainty than bright galaxies, their intrinsic redshift 
distribution may be wider. In this case, smoothing with 
the same Az will not sufficiently broaden the redshift dis- 
tribution of faint galaxies, resulting in an underestimate 
of the correlation length of faint galaxies for some ob- 
served value of A-u,. This will introduce an artificial trend 
of increasing clustering with increasing brightness. We 
note that many other studies of galaxy clustering which 
use the same redshift selection function for all galaxy 
samples-may be subject to this effect. 

A more appropriate redshift distribution for use in 
eq. 10 may be had by summing the redshift probability 
distribution P(z) (see §2) for each galaxy in the sample. 
To the extent that our Monte Carlo simulations provide 
an accurate estimate of P{z), this strategy circumvents 
the problem of choosing a smoothing width. Addition- 
ally, we compute ^ P{z) separately for each sub-sample 
under consideration, thereby reducing the problem of dif- 
ferential redshift uncertainties. 

The smoothed 2;phot distributions are narrower than 
the (more realistic) ^ P{z) distributions that are used 
throughout this paper. We note that using these nar- 
rower distributions would reduce our estimates of ro by 
~ 15 — 25%, where the more clustered populations dis- 
play the smaller differences. In the case of DRGs, about 
half of the ^ 15% difference between these two estimates 
of ro comes from the z < 1.8 tail of the ^ P{z) distribu- 
tion, and the other half comes from the broader overall 
distribution at z > 1.8 (Fig. 1). 

If redshift interlopers are assumed to be randomly dis- 
tributed they will dilute observed angular correlation by 
a factor (1 — /c)^, where fc is the contamination fraction. 
One method to account for interlopers is to estimate fc 
using the redshift probability distributions, and to calcu- 
late a contaminant-corrected ro by integrating over the 
range 1.8 < z < 3.5 in eq. 10. As the assumption of 
a random distribution is probably unrealistic, we choose 
instead to account for interloper's by integrating over the 
entire redshift range. The correlation lengths are similar 
regardless of the method used. 

3.3.2. Field to Field Variance 

The Landy & Szalay estimator has been shown to have 
approximately Poisson variance in the limit of zero clus- 
tering (Landy & Szalay 1993), but a clustered population 
is expected to show covariance between the radial bins 



(Bernstein 1994). The deep NIR MUSYC survey consists 
of only two independent fields (HDFSl/2 and 1030), so 
field to field variations will present an additional source 
of error. Moreover, uncertainty in the integral constraint 
is not correctly accounted for by Poisson statistics. Here 
we estimate the confidence intervals of our results with 
simulated data sets. Our approach is similar in spirit to 
that described by Daddi et al. (2003). 

We construct a clustered population using outputs 
from the public GallCS simulations (Hatton et al. 2003). 
GallCS uses cosmological N-body simulations to trace 
the growth and merging of dark matter halos, and a 
semi-analytic approach to follow the formation and evo- 
lution of galaxies within the halos. The simulations use 
8.272 X IO^Mq particles in a 100/i~^Mpc simulation box, 
Qm = 0.333, Oa = 0.667, h = 0.667, and as = 0.88. The 
GallCS outputs are available in convenient 'observing 
cone' catalogs (Blaizot et al. 2005) that mimic what an 
observer would see in a simulated universe. The limited 
size of the simulation box requires replication of galax- 
ies within the observing cones. Although the observing 
cone geometry has been tuned to reduce replication ef- 
fects, precise measurements of galaxy clustering and the 
cosmic variance are hampered by replication effects. 

We construct 27 mock data sets from the eight Ideg^ 
GallCS observing cones. Each of these mock data sets 
has the same geometry and field sizes as the deep NIR 
MUSYC survey, and we measure the clustering of the 
simulated galaxies using the same methods. The results 
are used to estimate the 68% confidence range of the 
amplitude of the angular correlation function. More 
detailed characterizations of the uncertainties would re- 
quire larger simulations. The confidence range is a func- 
tion of both intrinsic clustering, which is adjusted by 
selecting galaxies with different halo mass, and surface 
density, which is adjusted by randomly removing galax- 
ies. For most of the samples here, aA„/A,^, w 55% — 70% 
significantly larger than the Poisson values alone. This 
does not imply that our results are only significant at 
the < 2a level, as populations with small Ayj will have 
small aA„ , and so populations with little or no clustering 
rarely show strong clustering (see also Fig. 1 of Daddi et 
al. 2003). We will return to this issue in §3.4.2 for the 
specific case of DRG vs. LEG clustering. 

The errors due to field to field variations are reduced 
when comparing populations of galaxies drawn from the 
same fields, so we qiiote Poisson uncertainties except 
where noted. The estimated total uncertainty in the cor- 
relation length due to both Poisson errors and field to 
field variance is given in the last column of Table 2. We 
note that most studies of galaxy clustering at high red- 
shift do not fully account for the effects of field to field 
variance. 

3.4. Results 

3.4.1. Angular and spatial clustering of galaxies with 
2 < -Zphot < 3.5 

Fig. 3 shows the angular correlation function for K- 
selected galaxies with 2 < 2;phot < 3.5 and K < 21. The 
correlation function is roughly consistent with a power 
law down to ~ 1", the approximate resolution limit of 
our survey. There is a slight indication of an excess on 
smaller scales: the amplitude of the best-fitting power 
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Fig. 3. — The angular correlation function, with the integral 
constraint correction applied to the data points, for galaxies with 
2 < 2phot < 3.5 and K < 21. The upper x-axis shows the trans- 
verse comoving distance at the median redshift z = 2.6. The 
dashed hnes show the best fit power laws over 8 < 200" (upper) 
and 40" < 9 < 200" (lower). The lower fit is preferred in or- 
der to reduce the effects of halo substructure in w{6). The solid 
curve illustrates the shape of the angular correlation function that 
is expected for dark matter halos; the larger values of the galaxy 
correlation function at small separations suggest that some halos 
host multiple galaxies. 

law = 1.9 ± 0.3 when the fit is restricted to 9 < 200" 
and = 1.3 ± 0.4 over 40" < 9 < 200". Lee et al. 
(2006) and Ouchi et al. (2005) also note that the best- 
fitting power law changes for smaller angular intervals for 
their sample of z ~ 4 LBGs. Although our sample is not 
large enough to trace the detailed shape of ^(6*), it is pos- 
sible that this is evidence of the small-scale excess that 
is predicted by halo occupation distribution models (e.g. 
Wechsler et al. 2001; Zheng 2004). However, as empha- 
sized by Adelberger et al. (2005a), the observation that 
approximate power law behavior extends to such small 
scales may itself be interpreted as evidence that galaxies 
share halos. The solid line in Fig. 3 shows the expected 
shape of the angular correlation function of dark matter 
halos Wh{9), derived using N-body simulations (§3.3.2). 
Halo exclusion effects force Wh{9) to flatten on smaller 
scales. In the case of one galaxy per halo, w{9) should 
follow Wh{9). If galaxies have a higher probability than 
halos of having close neighbors, then it is likely that some 
fraction of these galaxy neighbors reside in the same halo. 
We return to this point in §4. 

We invert the angular correlation function w{9) to de- 
rive the spatial correlation length ro using eq. 10. Re- 
stricting the fit to the angular range 40" < 9 < 200", 
we find To = 6.0l5'^ft.~^Mpc (comoving). For com- 
parison, the correlation lengths of i? < 25.5 optically- 
selected z ~ 2 — 3 BX galaxies and LBGs is ~ Ah'^Mpc 
(Adelberger et al. 2005a; Lee et al. 2006). The lat- 
ter use a power law slope of the correlation function 
7 ~ 1.6, whereas we use 7 = 1.8. Assuming 7 ~ 1.6 
increases the correlation length of MUSYC galaxies to 
ro = 6.6l} °ft."^Mpc. The larger correlation lengths 
for ii'-selected sample might have been expected, as K- 



bright galaxies have been shown to cluster more strongly 
than optically-selected K-iaint galaxies at z ~ 2 (Daddi 
et al. 2004; Adelberger et al. 2005b). This dependence 
on selection filter may reflect underlying trends with K- 
magnitude, color, mass, or other parameters; these issues 
are discussed in the following subsections (see particu- 
larly §3.4.3). 

3.4.2. Clustering as a function of color 

Fig. 4 compares the angular correlation function of 
DRGs with J - K > 2.3 to that of non-DRGs with 
J — K < 2.3, in the redshift range 2 < Zphot < 3.5 and 
with K < 21. The DRGs cluster more strongly than the 
bluer galaxies at large scales. The angular correlation 
function of DRGs is roughly consistent with a power law 
down to 4", which corresponds to the virial radius of 
a ~ 1O^^M0 halo. For the blue galaxies, w{9) remains 
consistent with a power law to ~ 1". The two w{9) data 
points at 6* < 4" correspond to 1 and 8 observed blue 
galaxy pairs respectively, whereas the DRGs have and 
2 pairs at these small separations. Extrapolating the 
power law fit to these scales, we would expect only 1 and 
4 galaxy pairs for the DRG sample. While galaxy pairs 
at separations approaching 40" may share the same halo, 
it is not clear whether or not they will eventually merge. 
However galaxies with these significantly smaller separa- 
tions (4" corresponds to a projected distance of 36 proper 
kpc at z = 2.6) should be interacting strongly and could 
be mergers in progress. The lack of close DRG pairs 
may therefore indicate that DRGs are undergoing few 
mergers. Whether this is because DRGs are the prod- 
ucts of recent mergers, or because merger-induced star 
formation makes the red galaxies bluer, or whether there 
is some other explanation, is unknown. Evaluating the 
significance of the lack of close DRG pairs is severely 
complicated by galaxy deblending issues in our ~ 0.9" 
FWHM images, and we will not discuss this issue fur- 
ther. We find ro = llAtHh-^Mpc for DRGs. Note 
that, if the full field to field variance is taken into ac- 
count, we estimate ro = ll.lt42/i~^Mpc. If the fit is 

performed over 9 < 200", then ro = U.Oti^gh-'^M-pc 
(Poisson errors only). These values are consistent with 
the ro = ISAtf^h-^Mpc given by Grazian et al. (2006) 

and ro = 14.5;3;7/i-iMpc given by Daddi et al. (2003), 
although the latter authors do not apply a photometric 
redshift cut. 

Fig. 5 shows the comoving correlation length as a func- 
tion of minimum J — K color threshold. We confirm the 
previous result of Daddi et al. (2003) that the redder 
galaxies cluster more strongly, even though their result 
was derived using a single 4.5 arcmin^ field, and they 
measure w{9) at 9 < 70" (see their Fig. 8). It should 
also be noted that their sample reaches 2 — 3 magnitudes 
deeper than ours, and it is not obvious that the same 
trends should hold over such a wide luminosity range. 
There is a slight trend of increasing median redshift with 
increasing J — K color, but the difference is « 0.1 over 
the range of colors studied here, so it is unlikely that the 
relationship between color and ro is solely due to red- 
shift evolution. We note that alternate galaxy colors, 
such as i? — if, are also strongly correlated with cluster- 
ing (Fig. 6). LBGs and BX galaxies have a correlation 
length ro w 4/i^^Mpc (Adelberger et al. 2005a; Lee et 
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Fig. 4. — The angular correlation function for 2 < Zphot < 3.5 
galaxies that meet the J—K > 2.3 threshold for distant red galaxies 
(DRGs), and for non-DRGs. A small horizontal offset has been 
applied to the blue points for display purposes only. The power 
law fits to w{9) are performed at 40" < 9 < 200" to reduce the 
effects of halo substructure in w{9). 

al. 2006), lower than the value for the bluest threshold 
shown in Fig. 5, although the brightest i? < 24 LBGs 
reach tq = 7.8 ± O-Sft-^Mpc (Lee et al. 2006). The me- 
dian J—K color of z ~ 3 LBGs is ~ 1.6 (Shapley et al. 
2001), and very few LBGs/BX galaxies reach the reddest 
thresholds (Reddy et al. 2005; van Dokkum et al. 2006). 

We have established the significance of the increased 
clustering with color in several ways. Splitting our 
sample at the median color, J — K ^ 2.17, we find 



Fig. 5. — The comoving correlation length for 2 < z^^^i < 3.5 
galaxies redder than the J — K color threshold. 



^Mpc and = G.ltj.s^'^Mpc for the 



ro = ll.{)t\lh- 
redder and bluer sample, respectively. We then randomly 
split the sample of iC-selected galaxies in two repeatedly, 
measuring the correlation length each time. The corre- 
lation length reaches as high as rp = ll/i~^Mpc only 
~ 4% of the time, indicating that we have established 
the stronger clustering for redder samples at the ~ 96% 
level. We have also used the simulations described in 
§3.3.2 to see how often a population with the same cor- 
relation length as LBGs, but number density and redshift 
distribution similar to what we infer for DRGs, can have 
a measured correlation length as high as that observed 
for DRGs as a result of field to field variations; we found 
that this only happens ~ 5% of the time. Futhermore, 
we have verified that the increase in clustering with color 
is not driven by any one of our three ISPI fields by re- 
peating the clustering measurements three times, each 
time removing one of the fields; although the exact val- 
ues of the correlation length vary, the relationship be- 
tween clustering and color is always present. Finally, we 
recall that the 'total' uncertainties, which include the es- 
timated contribution from field to field variations and 
are presented in Table 2, overestimate the uncertainties 
when comparing correlation lengths of galaxies that are 
drawn from the same fields. 

The increasing clustering with color indicates that a 
color-density relationship was in place at z > 2. In the 
local universe, this relationship is understood as an ef- 




FlG. 6. — The comoving correlation length for 2 < Zphot < 3.5 
galaxies redder than the Rab — Kvega color threshold. 

feet of higher metallicity and higher stellar ages in the 
densest regions; both effects may play a role at high red- 
shift (Forster Schreiber et al. 2004; van Dokkum et al. 
2004; Shapley et al. 2004). The high dust obscuration 
associated with vigorous starbursts also contributes to 
the red colors of many K-selected galaxies (e.g. Forster 
Schreiber et al. 2004; Labbe et al. 2005; Webb et al. 
2006). It is entirely possible that the dusty and the "red 
and dead" galaxies (Kriek et al. 2006b) have different 
clustering properties, but the strong relationship between 
J—K and tq suggests that neither of these populations is 
weakly clustered. Disentangling the relationship between 
clustering and star formation for red galaxies would likely 
require large fields with Spitzer Space Telescope observa- 
tions (Labbe et al. 2005). 
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Fig. 7. — The upper plot shows the comoving correlation length 
for galaxies fainter than a A'-magnitude threshold. The grey points 
show the results from the HDFSl and 1030 fields only; these points 
have been offset in the horizontal direction for clarity. The lower 
plot shows that the fraction of red galaxies (defined as galaxies 
redder than the median, J — K > 2.17) does not change strongly 
over this magnitude range. 

3.4.3. Clustering as a function of apparent magnitude 

Fig. 7 shows the relationship between correlation 
length and minimum X-magnitude. There is a small 
trend with the fainter galaxies clustering more strongly, 
but the significance of this effect is low; removing either 
the HDFSl or HDFS2 fields from this analysis eliminates 
this relationship, while removing 1030 actually increases 
it. We conclude that the data do not suggest a strong 
relationship between K and tq. This contrasts with re- 
sults from z ~ 2.3 BX objects, which show clustering 
that increases strongly with K ( Adelberger et al. 2005b) ; 
we comment further on this below. The bottom panel of 
Fig. 7 shows that the fraction of red galaxies — where here 
we characterize a galaxy as red if it has J — K larger than 
the median — does not vary significantly with K . Com- 
bined with the result that rg correlates with J — K, it 
appears that color, and not if-magnitude, is the primary 
determinant of clustering strength. 

Next we split the sample into two populations using 
an apparent optical magnitude cut i? = 25, which is ap- 
proximately equal to the median total i?-magnitude, and 
is 0.5 magnitudes shallower than the limit used for z ^ 3 
LBGs (Steidel et al. 2003). Fig. 8 shows that the opti- 
cally brighter ii'-selected subsample clusters less strongly 
than the fainter subsample. This result is at odds with 
several studies of optically selected galaxies, including 
BX objects at z ~ 2.3 (Adelberger et al. 2005a)^, LBGs 
at z 3 (Giavalisco & Dickinson 2001; Foucaud et al. 
2003; Adelberger et al. 2005a; Lee et al. 2006), and LBGs 
at z - 4 and z ~ 5 (Lee et al. 2006; Ouchi et al. 2004), aU 
of which display stronger clustering with increased rest- 

^ It appears that for BX objects clustering strength increases 
with i?-brightness, ii'-brightness, and R — K color (Adelberger et 
al. 2005a, b); it follows that the faint blue BX objects are the least 
clustered. 



Fig. 8. — The comoving correlation length for ij-selected (blue 
symbols) and /C-selected galaxies (red circles), shown at represen- 
tative i?-magnitude values. The large triangles are for 2; ~ 3 LBGs 
from Lee et al. (2006). The small triangles are LBGs at 2 ~ 3, 
the squares are BX objects at z ~ 2.3, and the asterisks are BM 
objects at z ~ 1.7, all taken from Adelberger et al. (2005a). The 
LBGs and BX objects show increasing clustering with increasing 
i?-brightness, whereas the BM objects show no apparent trend, 
and the it-selected MUSYC galaxies appear to show the opposite 
trend. 

frame UV luminosity. It is interesting that our brighter 
subsample has a correlation length that agrees well with 
R-selected samples in the same magnitude range. Thus 
our results suggest that the X-selected galaxies that are 
below the limits of current i?-selected surveys are the 
most strongly clustered. 

We note that the median magnitudes of our two sub- 
samples are similar-with the optically- faint sample 0.1 
magnitudes brighter than the optically-bright sample- 
and the overall distributions of if-magnitude are also 
similar. So the anti-correlation between tq and R- 
brightness may simply be a manifestation of the correla- 
tion between tq and the R — K color (Fig. 6). Similarly, 
the observations of Shapley et al. (2005) suggest a re- 
lationship between R — K and K for BX objects, and 
Adelberger et al. (2005b) speculate that the observed re- 
lationship between ro and K for their sample may refiect 
an underlying correlation between tq and R — K. Thus 
the results for both X-selected galaxies and optically- 
selected galaxies suggest that color may be the most im- 
portant driver of clustering strength. If this is the case, 
the difference in colors between these two populations 
may explain the difference in clustering properties, as 
iiT-selected galaxies tend to be redder. 

3.4.4. Clustering as a function of stellar mass 

To investigate the relationship between stellar mass 
and clustering, we estimate the mass of MUSYC galax- 
ies by fitting Bruzual & Chariot (2003) models to the 
observed photometry at fixed Zphot- We assume a r = 
300Myr declining star formation history, solar metallic- 
ity, and a 0.1 — IOOMq Salpeter (1955) initial mass func- 
tion. 

Fig. 9 shows that the angular correlation function 



Clustering of i^T-selected Galaxies 



9 



of galaxies with > IO^^Mq and for galaxies with 
Af* < 10^^ M© are very similar. The top panel of Fig. 10 
shows the correlation length vs. minimum stellar mass 
threshold. There is no clear trend. The bottom panel 
of Fig. 10 shows that the fraction of red galaxies is an 
increasing function of mass threshold. It was shown 
above that rg increases strongly with J ~ K color; from 
this it might be expected that clustering would increase 
strongly with mass, because the most massive galaxies 
tend to be red. Moreover, the median masses of the R- 
faint and i?-bright samples in Fig. 8 are 6.2 x IO^^Mq 
and 1.6 x IO^^Mq, again suggesting a possible relation- 
ship between mass and rg. So it is interesting that we do 
not observe a clear relationship between clustering and 
stellar mass, although it must be noted that the data 
shown in Fig. 10 have large error bars. 
As discussed by van Dokkum et al. (2006), the majority 
65% in the current sample) of > lO^^M© galaxies 
at z > 2 are DRGs. With such significant overlap be- 
tween the massive and red galaxies, it is to be expected 
that they have similar correlation lengths. However, we 
find ro = 5.9t2;4^"^Mpc for Af, > 10"Mq galaxies, and 
ro = ll.lli[;^/i-iMpc for DRGs. This indicates that ei- 
ther the massive non-DRGs have a very low correlation 
length, or that the less massive DRGs have a very high 
correlation length. Our sample is not large enough to in- 
vestigate each of these sub-populations individually, but 
we do note that the correlation length for galaxies that 
are both massive and meet the J — K > 2.3 criterion for 
DRGs is ro = 7 .9^2^h~^Mpc, intermediate between the 
massive and DRG samples. This may be evidence that 
the low-mass DRGs are highly clustered and that the 
high-mass blue galaxies are less clustered. Relative to 
the median high mass DRG, the median low mass DRG 
is fainter in the NIR (K = 20.8 vs. 20.3), brighter in the 
optical {R — 25.8 vs. 27.1), but has a similar NIR color 
(J ~ K = 2.6 vs. 2.7). We have verified that it is not 
K-faint DRGs which contribute so strongly to the clus- 
tering but rather it is the low-mass DRGs by measuring 
ro = 10.5tlih~^Mpc for K < 20.7 DRGs. This result is 
analogous to conclusions from the local universe, where 
low-mass red galaxies and high-mass red galaxies inhabit 
the densest environments (Hogg et al. 2003; Kauffmann 
et al. 2004), and to the fact that, among massive galaxies, 
there is a strong relationship between correlation length 
and optical color (Li et al. 2006). Additionally, Kauff- 
mann et al. (2004) show that the stellar mass of galaxies 
is not a strong function of halo mass in the most massive 
halos. 

It should be noted that, while we are approximately 
complete for galaxies with Af* > 10^^ Mq (van Dokkum 
et al. 2006), we will be very incomplete for less mas- 
sive galaxies. Using stellar mass estimates from the ul- 
tradeep FIRES MS 1054-03 field (Forster Schreiber et 
al. 2006), we estimate ^ 65% completeness for galaxies 
with Af* > lO^^ '^AfQ. Scatter in our mass measurements 
may also obscure any relationship between ro and mass. 
Significantly deeper data are needed to study the depen- 
dence of rp on mass for a complete sample. 

4. RELATIONSHIP BETWEEN GALAXIES AND DARK 
MATTER HALOS 

In the halo model of galaxy formation, the galaxy cor- 
relation length is related to the mass of dark matter halos 
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Fig. 9. — Same as Fig. 4, for mass-selected samples. 
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Fig. 10. — The comoving correlation length as a function of stellar 
mass threshold. The lower plot shows the fraction of red galaxies 
(defined as galaxies redder than the median, J — K > 2.17) as a 
function of mass. 

(Mo & White 1996). In this section we constrain the halo 
masses and occupation numbers of the various subsam- 
ples of MUSYC galaxies using the measured correlation 
lengths and number densities. 

4.1. The number density and bias of dark matter halos 

To investigate the relationship between the MUSYC 
if-selected galaxies and dark matter halos, we use the 
halo mass function of Sheth & Tormen (1999) that is 
derived from fits to large N-body simulations 

dnu J 1 \ [2 p dv' f-v'^\ , , 
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where v' = ^JaSd cr{M , z), and the constants 5c ~ 1.69, 
a = 0.707, A w 0.322, q = 0.3, and p is the current mean 
mass density of the universe. We calculate the relative 
mass fluctuations in spheres that contain an average mass 
M as 

(T{M,z)=D{z)a{M,Q), (14) 

where D{z) is the growth factor for linear fluctuations 
given by Carroll ct al. (1992), and cr(Af, 0) is calculated 
using a scalc-frec n — 1 initial power spectrum and the 
transfer function of Bardccn ct al. (1986). 

The linear halo bias is calculated using the function of 
Sheth, Mo, & Tormen (2001) 



1 



j/'2c + 6(l-c)(l-c/2) 

(15) 

where & = 0.5 and c = 0.6. Further details can be found 
in e.g. Mo & White (2002). 

Several definitions of the bias-which relates the clus- 
tering of objects to that of the overall dark matter 
distribution- in terms of observable quantities appear in 
the literature. We choose 



(16) 



where <ys{z) is the variance in 8/i~^Mpc spheres and is 
calculated analogously to eq. 14. If the galaxy correlation 
function ^ (r) is a power-law of the form eq. 9 then it can 
be integrated to give the relative variance 



72 



(3 - 7)(4 - 7)(6 - 7)2T VS/i-^Mpc 



To 



(17) 



(Peebles 1980, §36, §59). 

We model the simple case of one galaxy per halo above 
a minimum halo mass threshold, i.e. a halo occupation 
number of 1. More detailed models, such as setting the 
halo occupation number equal to a power law above some 
mass threshold (e.g. Wechslcr ct al. 2001) are beyond the 
scope of this paper. As we measure the bias of the galaxy 
samples at scales larger than ^ IMpc, we can associate 
the observed bias with the linear bias of the host halos 
calculated with eq. 15, thereby providing an estimate of 
the halo mass. 

Fig. 11 shows the average halo bias (weighted by the 
number density) and number density as a function of halo 
mass threshold at the median redshift, z ~ 2.6. The Ict 
bias range for DRGs, as well as the minimum mass and 
the number density of halos with the same bias values, 
is illustrated by the shaded region. From this figure we 
can read off the mass threshold of halos that have the 
same bias as DRGs, Af;, « 1.5 - 3 x W^Mq. The larger 
sample of ii" < 21 galaxies has « 1.5 - 5 x IO^^Mq. 
For comparison, the LBGs and BX galaxies occupy halos 
with mass threshold Mh « IO^-^Mq and fa 
respectively (Adelberger et al. 2005a) . 

4.2. The halo occupation number 

In the local universe, halos with mass < lO^^A'/0 tend 
to have only one bright galaxy, whereas > 10^** M© halos 
may contain dozens (e.g. Kauffmann et al. 2004). There 
is also evidence of galaxies sharing halos at high redshift 
(Daddi et al. 2003; Zheng 2004; Adelberger et al. 2005a; 
Ouchi et al. 2005; Lee et al. 2006). We define the mean 



halo occupation number Noccup as the ratio of galaxy 
number density to the number density of host dark mat- 
ter halos. Occupation numbers greater than unity sug- 
gest that multiple galaxies can reside in a single halo. 

The simplest estimate for the number density of galax- 
ies comes from dividing the observed number of galaxies 
by the volume probed by our survey at 2 < z < 3.5. 
However there may be significant evolution of the ac- 
tual number density over this redshift range, and there 
is probable contamination by interlopers. We attempt 
to correct for these effects by estimating the fraction of 
the observed galaxies that lie at z ~ 2.6 using the red- 
shift distributions discussed in §2. We note that these 
two estimates agree to within the ficld-to-field variance 
within the survey, which is ~ 20%. The estimated 
number density of galaxies at z ~ 2.6 with K < 21 
is (5 ± 2.5) X lO^-^/i^oMpc"^. For DRGs, we estimate 
(2± 1) X lO-^ftfoMpc-^ GOODS-CDFS - the only other 
public field with size, depth, and multiwavelength cov- 
erage comparable to one of our fields - shows a lower 
density of massive and red galaxies than are present in 
the MUSYC survey, indicating that our field-to-field vari- 
ance may not be representative (van Dokkum et al. 2006; 
Grazian et al. 2006). Incompleteness and possible sys- 
tematics in photometric redshifts further complicate den- 
sity estimates; we therefore assign approximate 50% un- 
certainties. These number densities are consistent with 
estimates from the luminosity function (Marchesini et al. 
2006). 

Fig. 12 compares the number density and correlation 
lengths for various samples to the values that are ex- 
pected in the case of a one-to-one relationship between 
galaxies and dark matter halos. As noted by e.g. Adel- 
berger et al. (2005a), the observed properties of 2; ^ 3 
LBGs are roughly consistent with such a relationship. 
The same is true for the entire MUSYC i^-sclccted sam- 
ple, as well as the sample of A'/, > IO^^A'/q galaxies. 
However the redder galaxies deviate strongly from the 
expected relation, suggesting high occupation numbers. 
In particular, the number density of halos with the same 
bias as DRGs is 2.4 - 13 x IQ-^/ifoMpc"^, suggest- 
ing Noccup ~ 40^30- However, if the estimated field to 
field variance is taken into account (Table 2), the DRGs 
may have a correlation length as low as ro ~ 7/i~^Mpc, 
in which case Noccup ^ 1. The bright K < 20.5 BX 
galaxies from Adelberger et al. (2005b) also suggest high 
occupation numbers, although the very strong clustering 
in one of their observed fields may drive their result. 

4.3. Constraining the occupation number with the 
num,ber of close pairs 

In the previous subsection we estimated the occupation 
number Noccup using the standard method of comparing 
the observed number density of galaxies to the number 
density of their host halos. If many galaxies share the 
same halo, then the angular positions of galaxies should 
show strong 'clumps,' i.e. the correlation function should 
indicate large power on small scalers. Here; we obtain an 
independent estimate of Noccup by comparing the shape 
of the angular correlation function w{6) to the angular 
correlation function of dark matter halos Whifi). 

Although the angular and spatial correlation functions 
of dark matter halos are well-approximated by power 
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Fig. 11. — Upper: The number- weighted average linear bias 
as a function of minimum halo mass threshold at 2 = 2.6. The 
right axis shows the relationship between galaxy correlation length 
and the inferred large-scale galaxy bias. The shaded regions show 
the 1(T range of allowed tq, and the corresponding range in bias 
and halo mass. Lower: The comoving halo number density as a 
function of halo mass threshold. The shaded regions show the Icr 
range of mass and number density for the halos that host DRGs. 
The hatched region shows the observed number density of DRGs, 
illustrating that they are 40^3q times more numerous than the 
halos that host them. 
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Fig. 12. — The correlation length and the number density of dif- 
ferent populations of 2 > 2 galaxies. The LBG data are taken from 
Adelberger et al. (2005a). The ro for the _fC- bright BX galaxies is 
taken from Adelberger et al. {2005b). We estimate the number 
density of bright BX galaxies by using density of all BX galaxies 
given by Adelberger et al. (2005a) and applying a correction using 
the information given in Adelberger et al. (2005b). The solid line 
shows the approximate correlation length ro as a function of num- 
ber density that would be expected in the case of one galaxy per 
halo at 2 = 2.6; the upper and lower dashed lines show the same 
information at 2 = 3 and 2 = 2, respectively. Although the inferred 
halo mass for a galaxy population with a given ro will depend on 
the redshift of the galaxy population, the relationship between ro 
and number density for dark matter halos is not a sensitive function 
of redshift. 



laws on scales larger than ~ ryir , the probability of find- 
ing distinct halos with separation < r^r necessarily falls 
to 0; this sets ^h{r < r^r) ~ — 1- However, for a suffi- 
ciently large redshift selection window, 'Wh{0) will tend 
to remain flat on small angular scales because of projec- 
tion effects (see Fig. 3). The extent to which the galaxy 
correlation function w{0) differs from the Wh{0) on small 
scales can be used to constrain the occupation number 
(Wechsler et al. 2001; Bullock, Wechsler, & Somcrville 
2002; Adelberger et al. 2005a). 

The expected number of galaxies in the angular inter- 
val (0, 6„iax) around a randomly-chosen galaxy is 

^ = %^^L.(l + ^i'<e,_), (18) 

where Ug is the mean surface density of galaxies, and 
^<Sr,iax the value of the angular correlation function 
evaluated over the same angular interval (i.e. using eq. 1 
and a bin of width 0max)- This relation follows directly 
from the definition of the angular correlation function. If 
all galaxies are associated with halos, and if 9max is larger 
than the virial radius, the average number of galaxies 
within the same angular interval is the number of addi- 
tional galaxies within the host halo plus a contribution 
from neighboring halos, 

NKi f + NoccuprihTrel^^Jl + Wh,<e„„J, (19) 

where / denotes the average number of additional galax- 
ies in a halo that hosts at least one galaxy. Note that 
Noccup = ng/uh, so the right hand side of eq. 19 does 
not depend directly on Noccup- Combining eqs. 18 and 
19 gives an estimate of / that depends on both w{6) and 
Wh{0) over the interval {0,9max), as well as on Ug. Un- 
der the assumption that all halos above the minimum 
mass threshold host detectable galaxies, it is apparent 
that Noccup = if + !)■ If some fraction g of these halos 
do not host detectable galaxies, then Noccup is reduced 
by a factor (1 — g). There is some indication that g > 
for LBGs (Adelberger et al. 2005a; Lee et al. 2006). Here 
we make the simplifying assumption that g = 0, which is 
consistent with the Noccup ^ 1 measurements from the 
previous subsection, and note that the Noccup measure- 
ments of this section may be upper limits. 

We measure Wh,<i9^^^ directly from the observing cone 
output of the GallCS simulations (see §3.3.2), using halos 
with a large-scale angular correlation and redshift distri- 
bution similar to those inferred for our sample. Because 
the GallCS observing cone outputs do not specify the 
coordinates of halos, we instead use the coordinates of 
the most central galaxy in each halo when measuring the 
correlation functions. As shown in Fig. 3, Wh{9) flat- 
tens over the region 9 < 40". Our results are relatively 
insensitive to the details of this procedure, and the un- 
certainties are dominated by the uncertainties in w{9) 
and in the galaxy surface density rig, which we estimate 
from the observed variance among the MUSYC fields. 
We use the observed field-to-field variations to estimate 
the uncertainty in rig, and use 9max = 60" in eqs. 18 and 
19. 

We derive Noccup « 1.7 ± 0.3 for DRGs. Other galaxy 
subsamples have occupation numbers that are consistent 
with the DRG value to within Icr. In the previous sub- 
section we showed that the observed correlation lengths 
and number densities suggest Noccup 1 for the K < 2\ 
and ^ IO^^Mq galaxies, consistent with the values 



12 



Quadri et al. 



derived here. However, the value for DRGs is much less 
than the Noccup = 40I3Q that is inferred from the cor- 
relation length. This may indicate that our measured 
correlation length is an overestimate. Possible causes for 
this discrepancy are given in §6, however we note that 
these values are consistent if the estimated field-to-field 
variance in the correlation length is taken into account 
(§3.3.2 and Table 2). Previous studies have found corre- 
lation lengths even larger than ours (Daddi et al. 2003; 
Grazian et al. 2006), although their measurements may 
have been unduly influenced by small-scale structure in 
w{e) (§3.4, and Zheng 2004). 

In the previous subsection we showed that the number 
density and correlation length of K < 20.5 BX galaxies 
(Adelberger et al. 2005b) point toward very high occupa- 
tion numbers. However these authors indicate that their 
results do not change significantly if they include galaxy 
pairs at small separations in their analysis. This suggests 
that there is no evidence of a strong small-scale excess 
in their correlation functions, and therefore that bright 
BX galaxies also show a disagreement between their ob- 
served properties and the properties of dark matter halos, 
although the difference is not as significant as for DRGs. 

5. RELATING GALAXY POPULATIONS AT DIFFERENT 
REDSHIFTS 

5.1. Evolution of galaxy bias 

What are the z = descendants of the X-selected 
galaxies discussed in this paper? In the ACDM picture of 
structure formation, the large-scale distribution of galax- 
ies is determined primarily by the dark matter potential 
wells. Thus we can address the evolution of high-redshift 
galaxies by following the evolution of their host dark mat- 
ter halos. Fortunately, the dynamics of coUisionless dark 
matter particles are well-described by simple models or 
by cosmological N-body simulations. So while the com- 
plicated physics that dictates the evolution of the bary- 
onic component of galaxies (e.g. star formation, feedback, 
mergers) cannot be addressed by these models or simu- 
lations alone, we can constrain the bias and halo mass of 
the z = descendants. 

One way to investigate relationships between galaxy 
populations at different redshifts is to compare their ob- 
served bias. As the universe evolves with time, the dark 
matter becomes more clustered and the bias of a set of 
biased objects will decrease. The bias h{z) of a set of test 
particles evolved according to 



h{z) = l + {h{0)-\)/D{z) 



(20) 



(Fry 1996) where D{z) is the growth factor. It is im- 
portant to note that this equation does not account for 
merging and the evolution of the baryonic components of 
galaxies. Merging will play a role if galaxies in the dens- 
est (i.e. most biased) regions of space are more likely 
to merge than galaxies in less dense regions, thereby re- 
ducing the average bias of unique descendants. So it is 
possible that the bias will evolve faster than indicated by 
eq. 20. We refer to this as the 'galaxy-conserving' model 
of bias evolution. Fig. 13 shows tracks of bias evolution, 
along with the bias of different samples of galaxies (which 
we compute in a consistent way, using eq. 16). This fig- 
ure shows that the brightest LBGs at z ^ 3 (i? < 24; 
Lee et al. 2006) have a bias roughly consistent with the 
z ~ 2.6 K <21 galaxies studied here. At higher redshift. 
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Fig. 13. — The evolution of bias with redshift. The tracks show 
the evolution of bias calculated using the galaxy-conserving model, 
eq. 20. The filled red circles are based on this work, while the blue 
symbols are for optically-selected galaxies at various redshifts and 
black symbols are for local galaxies. The asterisks show the bias of 
BX galaxies from Adelberger et al. (2005a, b). The triangles show 
the bias of 2: ~ 3 LBGs from Lee ct al. (2006). The diamond and 
stars are for the 2 ~ 4 LBGs of Allen et al. (2005) and Ouchi 
et al. (2004), respectively. Open circles are from Zehavi et al. 
(2005). The square represents the richest cluster sample analyzed 
by Bahcall et al. (2003). In all cases, we have calculated the bias 
using eqs. 16 and 17. 

only the brightest z ^ A LBGs (i' < 24.8; Ouchi et al. 
2004; Allen et al. 2005) show biasing consistent with the 
lower-redshift DRGs, but the fainter z ^ 4 LBGs are not 
consistent. 

Fig. 13 shows that all of the high-redshift samples dis- 
cussed in this paper, including LBGs, evolve into highly- 
biased populations at z ~ 0. This point has been 
made previously for the optically-selected populations 
(e.g. Baugh et al. 1998; Ouchi et al. 2004; Adelberger 
et al. 2005a). Among the known galaxies at z ~ 2, it 
appears that only those that are faintest in K can be 
progenitors of typical field galaxies. 

5.2. Evolution of halo mass 

The preceding analysis only illustrates the bias of the 
descendants of high-redshift galaxies. If the galaxies 
within a given population follow different evolutionary 
paths, then the descendants will have diverse properties, 
and the average bias will have limited interpretive value. 
Knowledge of the range of environments or halo masses of 
the descendants is more meaningful. For a more detailed 
investigation of the z = descendants of z > 2 galaxies, 
we track the growth of dark matter halos using cosmolog- 
ical N-body simulations. We choose the GIF simulation 
for its size and mass resolution, and for the publicly avail- 
able halo catalogs and merger trees (Frenk et al. 2000). 
This simulation uses flm — 0.3, J^a — 0.7, h — 0.7, 
r = 0.21, and ag = 0.9. The linear size is Ul.Sh-^ co- 
moving Mpc. A minimally-resolved halo consists of 10 
particles with mass 1.4 x lO^°/i~^M0, but we note that 
the merging histories of halos less massive than ^ 100 
particles may be inaccurate (Kauffmann et al. 1999). 
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In these simulations, the position of a central galaxy 
within a halo is given by the position of the most-bound 
dark matter particle. When halos merge, the central 
galaxy of the most massive progenitor becomes the new 
central galaxy while the central galaxy of less massive 
progenitors, as well as any progenitor satellites, are kept 
as satellites in the descendant halo. We assume that the 
galaxies in the present sample begin as central galaxies 
in z ~ 2.6 halos more massive than the threshold masses 
given in the previous subsection, and follow the positions 
of these galaxies to z = 0. Occasionally satellite galaxies 
are ejected during mergers and may not be contained in 
any of the simulated halos at later times, however this 
effect occurs rarely for the massive halos considered here, 
and may be safely ignored. Additionally, we exclude ha- 
los near the edges of the simulation from analysis. 

We note that our treatment of halo evolution is differ- 
ent than that of many other authors, and these differ- 
ences may lead to contrasting conclusions. For instance, 
Grazian et al. (2006) employ the "merging model" of 
Matarrese et al. (1997) and Moscardini et al. (1998) as 
one method to study the clustering evolution of DRGs. 
In the context of this model, z = galaxies are consid- 
ered to be descendants of DRGs if they occupy halos that 
are more massive than the z > 2 DRG host halos. More 
halos will meet this mass threshold at z = than at 
z > 2, so many of these descendants enter the sample at 
intermediate redshifts, leading to a lower typical descen- 
dant halo mass. In this section we are concerned only 
with the direct descendants of z ^ 2.6 MUSYC galaxies. 

There is an essential ambiguity in interpreting the 
range of halo masses occupied by the descendants of high- 
redshift galaxies. The simulations show that-because of 
halo mergers-many z = halos host multiple descen- 
dants. However it is unclear whether the galaxies them- 
selves merge or whether they retain separate identities 
within a single halo. First we deal with the scenario 
in which the galaxies do not merge, as in the galaxy- 
conserving model discussed above. The red hatched re- 
gion in Fig. 14 shows the 68% range of host halo masses 
for descendants of high-redshift galaxies, as a function of 
halo mass at z ~ 2.6. The descendants of DRGs primar- 
ily occupy cluster-scale halos, with mass > 1O^'*M0. If 
our estimate of the correlation lengths for DRGs is cor- 
rect, then DRGs with if < 21 cannot be progenitors of 
the majority of local field early-type galaxies. It follows 
that DRGs exist in proto-cluster regions. The major- 
ity of LBGs also end up in group and cluster-scale halos. 
The black hatched region in Fig. 14 shows the mass range 
of halos occupied by descendants under the assumption 
that all galaxies within a single halo merge. In this case 
the most massive halos at z = 0, which host multiple 
descendants of high-redshift galaxies, are only counted 
once. The difference between the two hatched regions il- 
lustrates the importance of merging; for instance, > 50% 
of the halos that are inferred to host LBGs merge with a 
more massive halo between z ^ 3 and z ~ 0. This indi- 
cates either that LBGs tend to merge with more massive 
galaxies to form the brightest central galaxy in a halo, 
or that LBG descendants are satellites rather than the 
brightest central galaxy in these halos. It may be pos- 
sible that some of the progenitors of the low-mass red 
galaxies that exist in dense regions (e.g. Hogg et al. 2003) 
are LBGs, although this would require star formation to 
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Fig. 14. — The evolution of halo mass with redshift. The x- 
axis shows the halo mass at the redshift of observation, and the 
y-axis shows the mass of the descendant halos at z=0. The red 
hatched region indicates the predicted 68% halo mass range for 
the descendants of z ^ 2.6 galaxies under the assumption of no 
galaxy mergers. The black hatched region indicates the mass range 
under the assumption that all galaxies that share halos merge. See 
the text for details. The thick solid curves show the median halo 
masses under each of these scenarios. The horizontal lines at the 
bottom of the figure mark the range of halo masses inferred for 
LBGs (Adelberger et al. 2005a), K < Tl galaxies, and DRGs at 
the epoch of observation. 

cease shortly after the epoch of observation. 

6. SUMMARY AND DISCUSSION 

We have used SOOarcmin^ of UBVRIzJHK imag- 
ing from the MUSYC survey to study the angular and 
spatial correlation functions of if-selected galaxies with 
2 < Zphot < 3.5 and K < 21. The correlation 
length for this sample is vq = 6.0l5:?^~^Mpc, ~ 50% 
larger than optically-selected galaxies at similar redshifts 
(Adelberger et al. 2005a; Lee et al. 2006). The clustering 
of galaxies increases strongly with J—K and R—K color, 
and the J — K > 2.3 population of distant red galaxies 
has ro = 11.1^^ l/i'^Mpc. Our results for DRGs are 
lower than previous results (Daddi et al. 2003; Grazian 
et al. 2006). This may be partially due to our smaller 
uncertainties. Additionally, previous studies were only 
able to constrain the correlation function on small scales, 
where the signal may be strongly affected by galaxies that 
share the same halo (Zheng 2004; Ouchi et al. 2005; Lee 
et al. 2006). In contrast, we perform fits to the correla- 
tion function on scales larger than the typical halo virial 
radius where this effect is reduced. 

Nevertheless we confirm the basic trend indicated by 
previous studies, in which red galaxies are much more 
strongly clustered than optically-selected galaxies. More- 
over, this trend is not simply due to the fact that the 
red galaxies used in the previous studies were selected 
in the A'-band: clustering increases strongly with color 
even within our iiT-selected sample, while there is no 
significant relationship between i^'-magnitude and color 
or between X-magnitude and correlation length. These 
results suggest that a color-density relationship was in 
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place at z > 2, as is also indicated by the properties 
of optically-selected galaxies in different environments 
(Steidel et al. 2005). Whether this relationship is driven 
by galaxies that are red because of dust obscuration, or 
because of low specific star formation rates, remains to 
be seen. 

Cucciati ct al. (2006) and Cooper ct al. (2006) inves- 
tigate the fraction of red galaxies as a function of local 
galaxy density, finding that the color-density relation ex- 
tends only to 2;^ 1.3—1.5. This apparent contradiction 
with our results may be due to the difference in measure- 
ment techniques (fraction of red galaxies vs. correlation 
lengths). These studies may also be dominated by galax- 
ies that are less massive than the ^ ^ 1.5 descendants of 
MUSYC galaxies, so there may not be a contradiction 
if the evolution of the color-density relationship is mass- 
dependent (Cooper ct al. 2006). Finally, these studies arc 
based on spectroscopy of galaxies that are selected in the 
optical, leaving the possibility that they are incomplete 
for the reddest galaxies at these redshifts. 

A color-density relation at 2; > 2 has implications for 
current galaxy selection techniques. Adelberger et al. 
(2005b) show that, among an optically-selected sample 
of 2; ~ 2.2 BX objects, there is a strong correlation be- 
tween clustering strength and K-h&nd magnitude. Ad- 
ditionally, Shapley et al. (2004) show that these galaxies 
are also massive, metal rich, and have high star forma- 
tion rates. Taken together, these results indicate that the 
ii'-bright galaxies found by optical surveys show similar 
properties to those uncovered by NIR surveys (e.g. Franx 
et al. 2003; Daddi et al. 2003, 2004; van Dokkum et al. 
2004, 2006; Kriek et al. 2006a). It has been argued that 
the primary difference between optically-bright massive 
galaxies and optically-faint massive galaxies may be that 
the former are observed during a chance period of unob- 
scured star formation (e.g. Shapley et al. 2005). However, 
the observed relationship between clustering and J — K 
OT R — K color suggest that this is not the complete 
explanation. For instance, among our X-selected sam- 
ple, it is the R-faint galaxies (i.e. those that could not 
make it into the sample of Adelberger et al. 2005b) that 
cluster most strongly. Additionally, the very high cor- 
relation lengths measured by Adelberger et al. (2005b), 
ro ~ lO/i^^Mpc for galaxies with K ^ 20.5, may not 
be representative because one of their four fields shows 
anomalously high clustering; the average of the other 
three fields is ro = 5 ± l/i~-'^Mpc (see their Fig. 2). Ap- 
plying the same K < 20.5, R < 25.5 selection criteria to 
the MUSYC sample, we find ro = A.^tHhr^Mpc. We 
conclude that the difference between massive optically- 
selected galaxies and massive red JT-selected galaxies is 
probably not "transient" , but is instead related to a fun- 
damental difference in the host dark matter halos. 

It is interesting to consider whether ro is primarily re- 
lated to stellar mass, or some other property. We use stel- 
lar population models to estimate the stellar masses, and 
do not find a significant relationship between mass and 
ro. This is interesting, as there is a strong relationship 
between color and mass; together, these results indicate 
that either low-mass red galaxies cluster very strongly 
or that massive blue galaxies cluster less strongly. We 
find some evidence supporting both of these conclusions. 
We discuss our findings in light of the recent results 
from optically-selected samples, finding that there is ev- 



idence that redder colors are associated with larger cor- 
relation lengths for both samples. For optically-selected 
galaxies, the observed relationship between decreasing R- 
magnitude and increasing clustering strength for LBGs 
has been taken to suggest a positive relationship between 
halo mass and star formation rate at early times (e.g. Gi- 
avalisco & Dickinson 2001). On the other hand, stellar 
population synthesis modelling has not uncovered a re- 
lationship between rest-frame UV luminosity and stellar 
mass for typical optically-selected galax;ies (Shapley et 
al. 2001, 2005, but see Papovich et al. 2001 for a dis- 
cussion of a fainter sample). Similarly, Adelberger et al. 
(2005b) show that the clustering of optically-selected BX 
galaxies is strongly related to iiT-brightness. While there 
is a correlation between K and stellar mass (Shapley et 
al. 2005), there is also a relationship between K and star 
formation rate (Reddy et al. 2005; Erb et al. 2006). So 
while firm conclusions cannot be drawn from this evi- 
dence, it appears that clustering may not be determined 
by stellar mass alone for optically-selected samples. The 
evidence presented here for ii'-selected galaxies also in- 
dicates that color, rather than stellar mass, may be the 
primary determinant of rg. 

We compared the observed clustering length of 
MUSYC galaxies to that of dark matter halos, assuming 
the simple case of one galaxy per halo above a halo mass 
threshold. Galaxies with 2 < Zphot < 3.5 and K < 21 oc- 
cupy halos with M > 3 x IO^^Mq. The number density of 
galaxies is similar to the number density of halos, indicat- 
ing a mean halo occupation number of order 1. Galax- 
ies with redder colors reside in more massive halos, with 
DRGs residing in M > IO^^Mq halos. However DRGs 
are found to be more numerous than these halos, sug- 
gesting a mean halo occupation number Ngccup ~ 40^3o. 
If such large numbers of galaxies occupy the same halos, 
then it is expected that the angular correlation function 
w{9) will show very large values on scales corresponding 
approximately to the halo virial radius. We performed 
an independent estimate of the occupation numbers by 
comparing the small-scale values of w{0) to the values 
expected for dark matter halos; this yields occupation 
numbers ^ 1.5 — 2 for all samples studied in this paper, 
regardless of color and stellar mass. 

The cause of the discrepancy in the occupation num- 
ber of red galaxies is unclear. One possibility is that 
we have overestimated the correlation length of DRGs, 
although previous studies (Daddi et al. 2003; Grazian 
ct al. 2006) have found correlation lengths even larger 
than the value presented here. It is important to note 
that the high ro measured for DRGs is largely due to 
the integral constraint correction (§3.1); neglecting this 
correction decreases the power law amplitude of wlO) 
by ~ 50%, and decreases the correlation length from 
ro = ll.ltJ-j/i-^Mpc to ro = 7 .lt\lh-^MY>c. Larger 
fields are necessary to reduce the effect of the inte- 
gral constraint. We also note that occupation numbers 
greater than unity indicate that more complex halo oc- 
cupation distribution models in which the number of 
galaxies per halo depends on the halo mass-are necessary 
in order to accurately quantify the relationship between 
galaxies and halos (e.g. Zheng 2004; Lee et al. 2006). 
Regardless of the cause of this discrepancy, the relation- 
ship between color and clustering in our sample has been 
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established with ^ 96% significance. 

Finally, wc addressed the evolution of z > 2 galaxies. 
The descendants of our /sT-selected populations tend to 
occupy halos with masses 10^^ — IQ^^Mq, corresponding 
to the mass scales of groups. The reddest samples, in- 
cluding the DRGs, may occupy cluster-scale halos, with 
masses > IO^'^Mq. Even the descendants of the less clus- 
tered LBGs tend to reside in groups. It appears that 
only a small subset of the z > 2 galaxies that dominate 
current redshift surveys can be progenitors of typical X* 
field galaxies. 

An important caveat is that each galaxy 'population' 
will have rather heterogeneous properties, and it may 
be that discussing the evolutionary paths of population 
averages obscures important distinctions. For instance, 
Adelberger et al. (2005a) show that the correlation length 
of BX objects is ^ 4/i^^Mpc, and argue from this that 
their descendants should be elliptical galaxies. Adel- 
berger et al. (2005b) show that the BX objects with 
< 21 contribute most strongly to the clustering mea- 
surement; 40% of the BX objects are at K > 21.5, and 
have a correlation length ^ 2.5hr^Mpc. Presumably 
these fainter galaxies will evolve into a much less clus- 
tered population by z = 0. 

A principle limitation of the preceding analysis is the 
heavy reliance on photometric redshifts. Obtaining large 
numbers of spectroscopic redshifts for red J^-selected 
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galaxies has proven difficult on 6-lOm telescopes. While 
NIR spectroscopy yields a high success rate for determin- 
ing redshift for bright galaxies (Kriek et al. 2006b), the 
advent of multi-objects NIR spectroscopy will make the 
process more efficient. Another limitation of our study is 
the small galaxy sample. Recent clustering measurement 
of optically-selected galaxies are based on samples that 
are 1-2 orders of magnitude larger than the one presented 
here. The next generation of NIR detectors will enable 
the imaging of significantly larger fields to comparable 
depth, allowing for more precise measurements of galaxy 
clustering. 
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Clustering of ii'-selected Galaxies 



TABLE 1 

Galaxy correlation functions: fitting range 0" < < 200' 



Selection ^ 


Ngal 


A„(/3 = 0.8) 


ro(7= 1-8) 


r()(7 = 1.6) 


K <2l 


644 


1.9 ±0.3 


7 6+" " 

' -"-0.6 


■-'•-^-0.8 


J-K> 1.1 
J - K > lA 
J-K >1.7 
J-K >2.0 
J -K >2.3 


638 
614 
493 
381 
267 


1.8 ±0.3 
1.8 ±0.3 
2.3 ±0.4 

2.8 ±0.5 

4.9 ±0.7 


7 4+0.6 
'•^-0.7 
7 4+0.6 

o q+0.7 
».3_Q 7 

i2.oi?:^ 


8 9+''-'^ 
^■^-0.8 

O Q+0.8 

°-''-0.8 
10 6+°-* 

12 1+1° 

16.211:2 


R- K >2.9 
R- K>3A 
R-K>3.9 
R-K >4A 


626 
563 
444 
353 


1.9 ±0.3 
2.1 ±0.3 
2.9 ±0.4 
3.9 ±0.5 


7 t:+0.6 
' •■^-0.7 
o 1+0.6 
=•-^-0.7 
Q 6+0.7 
^•"-0.8 


10.2l°;| 
12.5l°| 
15.011:1 


R<25 
R>25 


341 
303 


1.9 ±0.5 
3.4 ±0.6 


7.3tl:« 

WAtli 


8.71- 
13.5l'lj 


K > 19.3 
K > 19.7 
K > 20.1 
K > 20.5 


620 
574 
480 
279 


2.1 ± 0.3 
2.4 ±0.3 

2.6 ±0.4 

2.7 ±0.7 


' ''-0.6 

8.6ig:? 

8.91^7 


9.6+?,-I 

io.4l°:« 
10.81°:^ 


log{M) > 10.4 
log{M) > 10.6 
log{M) > 10.8 
log{M) > 11.0 


616 
543 
429 
325 


2.1 ±0.3 

2.0 ±0.3 
1.8 ±0.4 

2.1 ±0.6 


o n+0.6 

7 (5+0.7 
' -"-0.7 

7.4«;« 


q 7+0.7 
'-0.8 

9.4i°:| 
9.311:^ 

9.611.^ 



''AH galaxies are selected using A" < 21 unless otherwise specified 



TABLE 2 

Galaxy correlation functions and bias: fitting range 40" < < 200' 



Selection * 




A„(/3 = 0.8) 


ro{7 = 1.8) 


ro(7 = 1.6) 


bias 


total variance'' 


X < 21 


644 


1.3 ±0.4 


6-ol?:? 


6-6Ii.i 


o ., + 0.5 
■^-•^-0.5 


+2.1 
-3.0 


J - ft: > 1.1 


638 


1.3 ±0.4 


"•^-1.0 


6.711:? 


o 4+0.5 
■^•^-0.5 


+2.0 
-3.0 


J - X > 1.4 


614 


1.4 ±0.4 


6 3+°-^ 
0-^-1. 


7.111:? 


q tr+0.5 


+2.1 
-3.0 


J-A" > 1.7 


493 


2.7 ±0.5 


q n+00 

='•"-'-0.9 


10.611:° 


A O+0.4 

4-8_o.4 


+2.6 
-3.4 


J-K >2.0 


381 


3.4 ±0.6 


10.2+1° 


i2.6ll:J 


^ q+0.5 


+2.9 
-3.8 


J-K >2.3 


267 


4.3 ±0.9 


n 1+1-3 

iJ^.i_1.4 


14-211:;! 


t- O+0.6 
tl.8_0 7 


+2.8 
-4.2 


R- K >2.9 


626 


1.5 ±0.4 


c c+0.9 


7.6+1-? 


■'-''-0.5 


+ 1.8 
-2.4 


R- K > 3A 


563 


2.3 ±0.4 


S K + 0.8 
°-f-0.9 


q q+"-» 


4 5+0-4 


+2.3 
-3.0 


R- K>3.9 


444 


3.2 ±0.5 


l0.2+'il 


12.4+1-1 


5 4+0.4 


+2.9 
-3.7 


R- K > iA 


353 


4.7 ±0.7 


12.511-;? 


15.5lli 




+ 1.6 
-4.3 


R<25 
R>25 


341 
303 


1.1 ±0.7 

3.2 ±0.8 


5.41^;^ 
lo.oil:^ 


5.9l-« 

11.911:? 


3.0+?-? 

^•^-0.7 


+2.3 
-5.4 
+3.0 
-4.2 


K > 19.3 
K > 19.7 
K > 20.1 
K > 20.5 


620 

574 
480 
279 


1.6 ±0.4 

1.7 ±0.4 
2.1 ±0.5 
3.3 ±0.9 


6.8i?:g 
7.oi?:° 

7.9+1-? 

lo.oii:^ 


7 7+1.0 
' ■ ' -1.1 

7.9+1-? 
9.0ll-i 

12-311:1 


7+0.4 
•^•'-0.5 

q O + 0.4 

5 2+°-'^ 


+2.1 
-3.0 
+2.1 
-3.0 
+2.6 
-3.4 
+3.1 
-4.4 


log{M) > 10.4 


616 


1.4 ±0.4 


6.3l?-? 


7 -,+1.0 
'-^-1.1 


o 4+0.5 


+3.1 
-4.4 


log{M) > 10.6 


543 


1.8 ±0.4 


7 9+0.9 
' -^-1.1 


8.2+1-? 


o g+0-5 
■^-='-0.5 


+2.2 
-3.1 


log{M) > 10.8 


429 


1.7 ±0.6 


7.itl;! 


8-2111 


O Q + 0.6 

cs.a_Q 7 


+2.5 
-3.5 


log{M) > 11.0 


325 


1.2 ±0.7 


5.91^:1 


6.2l|° 


Q q + 0.9 


+2.4 
-4.9 



''AH galaxies are selected using K < 21 unless otherwise specified 
''Estimated uncertainty in ro due to field to field variance. See §3.3.2 



